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Non-perturbative renormalization in kaon decays. 
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We discuss the application of the MPSTV non-perturbative method [l[ to the operators relevant to kaon decays. 
This enables us to reappraise the long-standing question of the AI — 1/2 rule, which involves power-divergent 
subtractions that cannot be evaluated in perturbation theory. We also study the mixing with dimension-six 
operators and discuss its implications to the chiral behaviour of the Bk parameter. 
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1. Introduction 

A quantitative theoretical understanding of the 
AI = 1/2 rule in K — > rnr decays has proven 
to be a formidable task since the calculation of 
hadronic matrix elements in the low-energy non- 
perturbative regime is needed. In the continuum, 
with an active charm quark and the GIM mech- 
anism at work, the operator basis given by 

= 7j[(sd) L (uu) L ± (su) L (ud) L ] -(u-> c).(l) 

In the framework of lattice QCD with Wilson-like 
fermions, the renormalization strategy is com- 
plicated by chiral symmetry breaking. In fact, 
the Wilson term induces the mixing of O ll 
with lower-dimensional operators, with power- 
divergent coefficients, which need to be sub- 
tracted non-perturbatively. Non-perturbative 
(NP) renormalization has witnessed a great 
progress recently jl] U . In this talk, we want 
to outline the strategy for the renormalization of 
the lattice operators Of L using the NP method 
(NPM) of 0. In the following, we shall concen- 
trate only on the octet component of Ofr , which 
we will denote with Oq In sec. g| we also 

reconsider the mixing with dimension-six opera- 
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tors [0J|] and its implications on the chiral be- 
haviour of Bk H , a problem also addressed in (J] 
using the Ward Identities (WI). We stress that 
the NPM of Q] and the WI are completely equiv- 
alent, provided the renormalization scale is large 
enough, since the renormalized operators with the 
good chiral properties are unique. Anyway, the 
NPM of yj must be used to calculate the overall 
costant, which the WI cannot determine. 

2. Renormalization strategy 

On the lattice, the general prescription to re- 
store the chiral symmetry broken by the Wilson 
term is to subtract all the operators of dimension 
equal or less than . The operators which need 
to be subtracted are dictated by the symmetries 
of the action: charge conjugation (C), parity (P) 
and s <-> d flavour switching symmmetry (S) [[10[ . 

There are different ways of calculating the K — * 
7T7T matrix elements, which correspond to different 
renormalization structures (2|(|. We considering 
here a general structure of the form 



o ± = z d 



z?o 3 



(2) 



being proportional to 
will be neglected in the following. 



where are the bare operators, Of, i = 
1, . . . , 4 are dimension-six operators of wrong chi- 
rality (cf. sec. ||), O5 is a dimension- five of the 
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form sB^yF^ v d (X^,, = or a^) and 3 is a 
dimension-three operator of the form sTd (T = I 
or 75). According to the NPM, the mixing Z's 
are determined by finding a set of projectors on 
the tree-level amputated Green functions (GF), 
with off-shell quark and gluon external states, the 
choice of which depends on the nature of the op- 
erators at hand. For the AI = 1/2 operators 
we choose the following set of external states: 
qq, qqg, qqqq, with the momenta given below in 
eq. (||). For each choice of external states, i.e. for 
each different set of GF, we need different type of 
projectors. Let us denote with P 3 the projector 
on the qq GF of the operator 3 , with P 5 the pro- 
jector on the qqg GF of the operator O5, and with 
Pj, j — 1, ... ,4 the set of mutually orthogonal 
projectors on the operators O,;, i = 1, . . . ,4 |||8|. 
Applying the projectors to the corresponding NP 
GF of the renormalized operators , with an 
appropriate choice of the external states, we re- 
quire that the renormalized operators be propor- 
tional to the bare operators, ± (fj,) oc O^(a) (up 
to terms of 0(a)), i.e. we impose the following 
renormalization conditions (trace over colour and 
spin is understood in the projection operation): 



P 3 {q(p)\d ± \q(p)) = 
PK{q{p-k)g(k)\d±\q(p)) 

q(p)q(p)) 



5\ 

■ ± 



P~(g(p)g(p)|0 ± 





= 0, j 



(3) 



where p and k denote the momentum of the ex- 
ternal quark and gluon legs. The system of equa- 
tions (^|) completely determines in a NP way the 
renormalization constants, as we have six condi- 
tions (non-homogeneous due to the matrix ele- 



ments of O , cf. eq. (|2j)) in six unknown mixing 
constants, Zf, i = 1, . . . ,4, Zf,Zf. 

Unfortunately, since solving eq. (|^) involves 
delicate cancellations between large contribu- 
tions, it may very likely result in a very noisy de- 
termination, even with large statistics. An equiv- 
alent strategy we can adopt is: 
1. We introduce an intermediate subtraction for 
the dimension-five and -six operators 



•Cf"0 S) » = 0,...,4, 



o 5 + c { 3 5) o 3 , 



(4) 



of - of 
o 5 

and determine the power-divergent mixing con- 



stants C3 and Cg 5 '' by imposing 



^3<?(p)|O i ± |?(p)) = 0, i = 0, 
P 3 (q(p)\0 5 \q(p)) =0. 



,4. 



(5) 



2. The finite mixing constants Z i and Z 5 , which 
in principle can be calculated in perturbation the- 
ory (PT), are then determined from the system 



P 5 (q(p - k)g(k)\ £\ ZfOf + Zf0 5 \q{P)) 
= P 5 {q(p-k)g(k)\d±\q(p)) 

Pj(q(p)q(p)\ zfof + z£d 5 \ q (p)q(p)) 



(6) 



Pj(q(p)9(p)\Of\q(p)q(p)), j = l, 
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3. Mixing with dimension-six operators 

In ref. the NP mixing with equal dimen- 
sional operators of the AS — 2 operator, which 
has the same Lorentz structure of Oq , has been 
studied, with the mixing limited to a basis of 
three operators, denoted here with Of, i = 
1,2,3, found by an explicit 1-loop calculation 
|ll"l| . A more careful analysis based on CPS0 
has shown that beyond 1-loop, the operators Oq 
are allowed to mix not only with O. , i — 1, 2, 3, 
present in 1-loop PT, but also with 

Of = {± ^ Nc 1) [(SS+PP-ItT)±(2 - 4)], (7) 

which is not present at the 1-loop level, where 

rr ee (^r^X^r^) and N c = 3. 



The operators 0\ 



0, . . . , 4 are Fierz eigen- 



states in Dirac-colour space with eigenvalue ±1, 
respectively || . We stress that Fierz transforma- 
tions are not a symmetry of the action, but only of 
the operators, thus they cannot determine the op- 
erator mixing. Nevertheless, since Fierz transfor- 
mations have been used in the perturbative cal- 
culations flllj , we have also used them to classify 
and reorganize the NP mixing. 

We have extended the method for the renor- 
malization of the AS = 2 operator presented in 
ref. [Q including the fourth operator O4 , and have 
determined the mixing constants Zi, i = 1, . . . , 4, 
and the overall constant Z AS=2 . We have per- 
formed the calculation using the same parameters 
of ref. 0, on a 16 3 x 32 lattice, at j3 = 6.0, with a 
hopping parameter n = 0.1425 for the SW-Clover 
quark propagator, in the lattice Landau gauge, 
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Figure 1. NP mixing Z's as a function of /x 2 a 2 . 
The solid (dashed) line is from SPT (BPT). 
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Figure 2. Chiral behaviour of B K . The NP Z's 
are taken at /x 2 a 2 = 0.96. 



but with higher statistics, i.e. on an ensemble of 
200 configurations. We have calculated the Z's at 
different renormalization scales /i 2 a 2 (cf. fig. Q). 
We note that Z% and Z4 are very well defined and 
almost scale-independent in a large "window" of 
/i 2 a 2 , whereas Z\ and Z 3 present a smaller win- 
dow, i.e. a more pronounced scale-dependence. 
Moreover, Z4 which is absent in 1-loop PT, is 
not negligible. We stress that the large fluctua- 
tions at small /i 2 a 2 do not spoil the validity of 
the NPM, since in that region the perturbative 
matching to a continuum scheme is not reliable, 
as for any NP lattice method 

Using the NP Z's of fig. |l|, the chiral behaviour 
of the Bk parameter studied in refs. can 
be revisited. Since the mixing with O4 starts at 
0(<7q) we do not expect drastic changes in the chi- 
ral behaviour. If the chiral behaviour were sensi- 
bly different, we would not trust the matching to 
the continuum which has an uneludable pertur- 
bative uncertainty. In fig. |^ the result using the 
NP Z's at /1 2 = .96 is shown. We will concentrate 
here on the novelties introduced by the NP mixing 
with the complete basis (N = 4) with respect to 
the basis found in 1-loop PT (N = 3) refer- 
ing to Q for a detailed analysis. Using the usual 
parametrization (O as=2 ) = a+/3m 2 K +^(p-q)+..., 



we have found that the use of the complete set of 
operators leaves the values of (3 and 7 almost un- 
changed compared to the N — 3 case, whereas 
the value of a is slightly lowered, passing from 
a N = 3 = 0.017(13) to a N=i = -0.004(15). We 
note that the two values of a are compatible with 
each other and with zero. Within the current 
statistical limitations forced by the thinning ap- 
proximation we can only state that the value 
of Bk, proportional to 7, is unaltered, while the 
correct chiral behaviour of the continuum, sig- 
naled by the vanishing of a, is recovered. These 
results are independent of the scale y, used, albeit 
not too small. We argue that this stability is due 
to the extremely clean determination of Z 2 ||] . 
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